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We  are  going  to  preaenit  here  a  method  of  eomputlng  the 
aerodynamle  eharaeterlstloa  of  a  .partidular  type  of  wing  of 
finite  span  under  the  oonditlnns  of  a.  streamlined,  gaseous 

flow.  ^  ^ 

If  wo  disregard  the  problem  bf  the  Interference  bet¬ 
ween  the  wing  and  the  fuselaaej  these  wings  can 
to  within  any  degree  of  precision  by  using  the  method  here  of¬ 
fered. 

No  examples  of  the  computation  of  such  wings  of  finite 
span  can  be  found  In  the  contemporary  literature 'on  gas  dyna¬ 
mics. 

1,  Shapes  of  Wings  Examined  ■ 

Let  us  suppose  that  the  airfoil  ABODA  In  the  drawing 
Inserted  represents  the  cross  . 

section  of  the  wing  at  the 
points  of  contact  with  the  fu-  _ 
selage,  while  one  airfoil,  cor-  ^ 

responds  to  the  profile  of  the 
same  wlzig  In  Its  oonsollo  cross 

section.  ^  — -Vt 

We  shall  examine  such  _ 

wings  AOHP,  whose  upper  and 
lower  sixrfaoes  combine  to  form 

the  side  surface  of  a  trimca-  JZ'  x 

tod  cone  with  the  bases  ABODA  ^  ' . .■  .  '  • 

and  FSHEF  and  the  apex  at  £. 

And  we  shall  assume  that  the  Fig.  1. 
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ioonsollo  crosa  section  waa  chosen  In  a  manner  which  either  ^ 
eliminated,  or  minimized  to  within  a  neghlglble  margin,  the 
effect  of  the  conaollo; cross  section.  Specifically,  this 
will  always  be  true  If  the  oonsollo  cross  section  degenera¬ 
tes  into  the  geometric  apex  JS  of  tbe  wing. 

2.  Differential  Equations  ai^  Boundary  Conditions. 


The  particular  shape  of ‘the  wings  ex^lned  here  Imme¬ 
diately  suggests  that  the  problem  of  streamlined  flow  should 
be  oonsldered  In  a  manner  reflecting  the  special  wj^-s^pe. 

For,  certainly,  all  points  of  any  straight  line  from 
the  apex  E  of  this  cone,  whose  side  surface  Is  the  surface 
of  the’  wliaig,  are  going  to  be  In  Identical  position  relative 
to  the  wlng^s  surxaoe  and  the  oncoming  flow,  ^  .  j.  ^ 

ItiS  natural,  therefore,  to  expect  that  under  stresm- 
llnlng  the  hydrodynamic  elements  shall  be  constant  along 

each  such  line.  ,  ^ 

Furthermore,  it  Is  natural  to  expect  that  the  resis¬ 
ting  surfaces  of  strong  disruption  of  the  stream  shall  also 
be  oonicftl  surfaces  with  apexi at  the  same  point  E  of  the 
streamllned-fluxed  cone  of  the  wing. 

To  verify  the  validity  of  these  expectations,  we  shall 
pass  from  the  rectangular  coordinates  x,  y,  z  to  different 
coordinates.  The  origin  of  that  system  will  coincide  vlth  E, 
but  we  shall  Introduce  two  new  coordinates  ^  and  according 
to  formulas 

^-T.  ’1-4*  (2.1) 


Let  us  now  consider  the  system  of  differential  eqi^— 
tions  of  motion,  and  transfonn  to  the  coordinates  |  and  e 
with  the  assumption  that  the  gas-dynamic  elements  depend  ex 
cluslvely  on  these  coordinates.  , 

Thus,  wo  obtain  the  following  system  of  differential 

equations: 

+  - IJ.  (2,2) 

iVj,  -  Cjl)  If  P  +  (Vy—Ca^^lgp (2*®) 
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Since  our  system  contains  none  of  the  z,  y,  s  ooordi<j> 
nateSf  but  only  their  combinations  In  terms  of  the  magnitu¬ 
des  %  and  n  ,  ve  can  sebh  solutions  depending  only  on  k  and 
n*  Tbls  is  sufficient  to  justify  our  expectations,  since 
it  la  obvious  that  with  the  actual  opnlclty  of  the  streamli¬ 
ned  surfaces  and  the  assumed  cohlolty  of  the  surfaces  of 
strong  disruption  in  the  gaseous  flux  we  ehall  be  able  to 
express  all  boundary  oondltlon.s  in  the  terms  of  these  new 
coordinates  |  and  I7  •  *  ,  ^ 

Let  the  equation  oft the  streamlined  cone  bet 


On  computing  the  direction  cosines  of  the  normal  to 
the  surface  of  the  cone  with  coordinate  axes  K,  y,  2,  and  on 
equating  the  normal  component  of  the  wing's  velocity  to  zero, 
we  shall  obtain  the  streamllnlxjg  condition  in  the  formt 

V,  -  -  (Vjt  —  iv,)  .  (2,8) 

The  condition  (2,8) ,  loj  expressed  in  terms  of  5  andy, 
which  justifies  our  expeotatibns. 

We  shall  pass  now  to  the  conditions  on  the  surface  of 
the  strong  disruption*  Let 


be  the  equation  of  this  surface. 

The  X  axis  will  be  directed  in  the  direction  of  the 
flow.  We  shall  compute  the  direction  cosines  of  -the  normal 
to  that  surface  and  substitute  these  in  the  expressions  hol¬ 
ding  for  disruption  surfaces*  The  relation  on  the  surface  of 
strong  disruption  will  have  the  forms: 

9^^  V,  -rl-i  4  f  . -v  .  ^1. - (2.»0> 


"  nn  v;  y  jt  * + wx'  -  n)*  “  ***]  ’ 

A  •  171^  [’'*’  r+  jr^ + ax'  -  n)’ ""  ^  •**]  * 
rn’*‘rnV^/ — 


(2.13> 

(2.I4> 


♦The  sign  in  formula  (2,12)  will  depend  on  the  choice  of  the 
direction  of  the  axis* 
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r“  The  Index  2  In  these  relations  corresponds  to  the 

drodynaalo  elements  after  the  disruption,  and  the  Index  1  to 
the  values  prior  to  thh  dlsx*uptlon« 

These  relations  (2,10)  -  (2,14)  are  again  expressed  In 
the  terms  |  and  t1  ;  consequently,,.  o\ir  expectations  are  fully 
Justified,  •  ,  .  ,  ,  ^  ... 

In  reference  to  all  the  above  facta,  conclude  that 
o\ir  Inveatlgatlons  must  be  applied  In  the  plane  of  these  two 
coordinates  (fe,-^).  In  this  plane  we  shall  have  to  ^tegrate 
the  system  of  equations  (2,2)  to  (2,6)  with  the  boundary  con¬ 
dition  (2,8)  on  the  profile  and  with  the  boundary  oondltions 
(2,10)  to  (2,14)  on  the  still  unknown  line  of  disruption. 

Should  bhe  system  of  the  equations  (2,2)  to  (2,6)  turn 
out  to  be  hyperbolic,  then  the  problem  can  be  solved  by  ap¬ 
plying  the  characteristic  method  In  the  numerical  integration 
of  l^perbollc  systems, 

3.  The  Type  of  the  System  of  Equations  and  Their  Charac¬ 
ter  let  Ics 


If  the  equation 


(3.1) 


is  the  characteristic  equation  of  the  system  (2,2)  -  (2,6), 
the  usual  methods  will  lead  to  the  following  equation  for  the 
determination  of  the  characteristic  directions  (characterized 
by  angular  coefficient 


A  «,  0,  -15"  ® 

0.  A,  0,  0 

0,  0.  A  0 

g,  I.  Ig-iv  0. 

0,  0,  0,  A,  — 


(3.2) 


where 


4  «  -  5(0,  -  ^>*)  +  (9,  -  yfOg) 


and  a  la  the  velocity  of  the  sound. 

Equation  (3,2)  Is  of  the  fifth  degree  In 

solving  It  we  shall  obtain  five  roots.  ,  7 

On  expanding  the  determinant  (3,2)  and  solving  the 
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raultlw?  oauatlon,  we  find  that  It  has  only  three  distinct 
root^beoause  one  of  these  roots  occurs  three  times. 

For  th©  roots  of  the  ohsraoterlstlc  deterailnftiitf  we 
obtain  the  expresslonb: 


(gj  "T 


r  "**  ***  ,1 

?El±I!I±jE.* 


:  e  V  le,  -  e4>»  -r  <».  - -f  ^ ± JOL .  (SJU 

(S) 

The  roots  of  the  character  let  lo  determinant  (^K-cr 
are  always  real,  while  the  roots  of  be  real  if 

the  inequality 

<2g  rig*y-±<»ryv>* + .<*>1  ~  >  a\  (3,5) 

holds.  In  that  case,  the  feystem  (2,2)  -  (2,6)  will  be  per¬ 
fectly  hyperbolic  and  the  characteristic  method  can  be  ap— 

obtain  an  idea  as  to  the  possible  realons  of  com¬ 
plete  hyperbollclty  for  the  system  (2,2)  -  (2,6),  we  shall 
consider  a  flow  hawing  constant  hydrodynamic  elements  and 
mowing  along  the  x  axis  with  the  supersonic  weloclty 

while  the  weloclty  of  sound  l8  «t.  ' 

In  this  case,  the  inequality  (3,5)  shall  be  trans¬ 
formed  into  the  lnequality_ 

(3.6) 


—  V  <  1. 
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This  inequality  hbiSs  in  the  region  between  the  two 
branches  of  the  hyperbola.*  The  region  of 
llclty  will  expand  as  Mach's  number  increases.  When  Mach 
number  is  less  than  1,  this  region  wanlshes. 

4.  Relation  on  the  Oharactorlstlos  and  Application  of 
the  Ghacacterlstio  Method 

The  relations  on  the  characterlstles  are  obtained  in 
the  usual  way.  The  roots  of  the  characteristic  equation  are 
substituted  sequentially  into  the  characteristic  system, 
oontalnim  partial  derlwatlwes  of  the  unknown  ^notions  with 
respect  to  one  independent  warlable.  The  resulting  relations 
do  not  determine  all  partial  derlwatlwes,  but  the  latter  can 
be  disregarded.  As  a  i*e8ult  of  the  exclusion  of  these  psJ^~ 
tial  derlwatlwes,  we  obtain  the  relations  on  the  characteris¬ 
tics. 

There  will  be  flwe  Independent  illations. 

We  will  get  one  relation  by  excluding  the  partial  de- 
plyatiwes  from  the  fozmiulas  which  we  obtained  by  substituting 
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th9  roots  in  the  oheraoterlftle  system:  and  three  for-] 

tne  tihe  roots^^),.^,*-  In  the  same  system. 


c 


;VuO  JTwWWO  %  (H  S.  r  -e**  ^ ^ T  i 

mulas  hy  substituting  the  roote{i^)x<^/»‘'  ^  r  . 

As  a  result  of  the  computations  we  get  the  follovimg 

relations  on  the  pharacteristiost 

- ijS - ••(*  +  V»«)(»x-“ H^»)l  +  «*«tj f.  (4.0 

(ejr**“tsi^<iejj4' (Wy  —  ■■  (4,2> 


(i|e,~le,)<le,  +  (vjr- 


(4.3> 

(4.4> 


=  The  relation  (4.1)  holds  along  the  characteristics 

with  the  Incllnatlnni?-^)/,^.  Foraulas  (4.2)  to  (4»4) 
hold  on  characteristics  with  the  Inclination 

Gtlven  the  five  relations  on  the  oharaoterlstlos  and 
the  boundary  conditions,  we  can  replace  these  by  formulae 
In  finite  differences  and  thus  find  the  five  unknown  func- 

Vv-e  Vvf  ^  fluid  pe 

The  order  in  which  these  functions  will  be  deter¬ 
mined  will  correspond,  In  genei*al,  to  the  order  In  which 
solutions  of  the  planar  whirlwind  problem  of  gas  dynamics 
are  found.  Accordingly,  see  any  adequate  textbook  on  gas 

d]7’31&H^O&m  ,  m 

No  difficulties  in  principle  will  be  encountered; 
thus,  there  Is  no  need  to  describe  in  detail  the  applica¬ 
tion  of  these  operations  to  a  practical  problem. 

5,  '/fings  of  Low  Conlolty  and  a  Slight  Arrow-Porm% 

If  the  wing  in  question  has  a  low  ponloity,  and  is 
oxay  slightly  arrow-shaped,  the  values  of  S 
velocity  regions  affectlbg  the  values  of  the  hydrod^amle 
elements  upon  the  surface  of  the  wing  shall  be  small. 

Let  us  assume  that  In  these  regions  the  values  of  % 
and  >7  are  such  that  their  squares  can  be  disregarded  com¬ 
pared  to  unity;  let  us  consider  otir  problem  under  this 

condlton.  . 

It  Is  clear  that  for  such  wings  in  these  regions 
where  the  values  of  the  hydrodynamic  elements  depend  on  ine 
wing,  the  order  of  the  magnitude  vg  and  of  Its  partial  de¬ 
rivatives  will  be  the  same  as  the  order  of  the  mwnltudes 
iL/t  /  and/t'y^/,  where  V;  Is  the  velocity  of  unperturbed  flow. 
'  Because  of  this  condlton,  we  can  consider  the  magni¬ 

tude  vx-^eaa  equal  to  Vx*  to  wlthto  the  ac^raoy  assumed  in 
the  discussion.  Simultaneously,^^  and^jJ^  win  1^®  ne¬ 
gligible  In  comparison  with  magnitudes  of  the  order  of  V;, 

°  ^  jftiriwp]  t.be  above  Into  consideration,  we  shall  get 

Wand  the  ma^ltude  |^“^»t  as  equal  to>^y 
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Finatead  of  the. system  (2,2)  -"(2,6),  the  system: 


.  fj. 

(6.IK 

(5.2) 

«jr ^  f  (^5  + ’i 

(S^V 

j  Wjr  J  ^  P  4- ^  -  0, 

(5.4> 

(5.5^ 

.  Equations  (5,1),  (5,2),  (5,4),  and  (5,5),  considered 
as  a  whole,  are  the  same  as  a  system  of  equations  for  the 
planar  motion  of  a  gas.  And  as  soon  as  these  equations  are 
Integrated,  the  Integration  of  (5,3)  can  be  performed  quite 
easily. 

The  streamlining  condition  (2,8)  Is  the  same  In  this 
ease  as  the  streamlining  condition  In  a  planar  flow  problem, 
axid  has  the  form 

. .  (6.«) 

Anif,  flxuCLly,  the  boundary  conditions  (2,10)  to  (2,14) 
assume,  under  the  assumptions  made,' the  followl^  forms: 


I  ' 

\v».7 

\V,.7 

i 

1 


P7 

h 


/•»  V 


W> 

.  (5.®) 

(5^ 
(5.10) 
(5.ll> 


The  conditions  (5,7)„  (5,8),  (5,10)  and  (5,11)  are  the 
same  as  the  corresponding  boundary  conditions  of  a  planar 
problem. 

Condition  (5,9)  must  be  used  to  determine  tbe  quantity 
Vw*  W®  obtain  the  quantity  by  Integrating  the  equation 
(5,3),  using  the  boundary  condition  obtained  from  condition 
( 5, 9  /  •  -  , 

Thus,  we  conclude  that.  In  order  to  solve  the  problem 

of  streamlining  for  a  wing  of  low  oonlcity  and  a  slight  arrow- 
form,  It  suffices  to  solve  the  problem  of  streamlining  of  a 
wing  of  an  infinite  spaa  having  a  profile  ^  s 
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Then,  we  transform  from  the  variables  ^  and  ft 
to  the  variables  x,  y,  a.  Thus,  the  computation  of  the  ae¬ 
rodynamical  character iatics  of  the  wings  under  consideration 
differs  only  negligibly  from  the  computation  of  the  aerodyna¬ 
mic  characteristics  for  wings  of  Infinite  span. 


E  N  D 
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